A complexity measure for continuous time quantum algorithms 
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We consider unitary dynamical evolutions on n qubits caused by time dependent pair-interaction 
Hamiltonians and show that the running time of a parallelized two-qubit gate network simulating the 
evolution is given by the time integral over the chromatic index of the interaction graph. This defines 
a complexity measure of continuous and discrete quantum algorithms which are in exact one-to-one 
correspondence. Furthermore we prove a lower bound on the growth of large-scale entanglement 
depending on the chromatic index. 
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I. INTRODUCTION 

At the moment, the most popular model of a quantum 
computer consists of the 2" dimensional Hilbert space 
Hn '■= (C^)®" of 'n qubits' as its memory space, and 
some one- and two-qubit gates as its set of basic transfor- 
mations (see e.g. [y ). There are several reasons for taking 
one- and two-qubit gates as the basic ones: Firstly, from 
the pure mathematical point of view, it is quite natural 
to look for a subset of the Lie group of unitary transfor- 
mations on the Hilbert space ?i„ generating the whole 
group. Obviously, one-qubit operations do not gener- 
ate the whole Lie group, the set of two-qubit gates does 
0. Hence there is no reason for taking more compli- 
cated transformations like three-qubit unitary operators 
as basic ones. Secondly, the model of two-qubit gates 
might be considered as the attempt to develop quantum 
computation in strong analogy to the theory of classi- 
cal devices: Building complex logical networks from two- 
bit gates is a successful concept of classical computation. 
The third reason stems from physics. Unfortunately it is 
only 'a little bit' convincing: From the fundamental point 
of view, n particles interact always in the form of pair- 
interactions, i.e., the total Hamiltonian H of the system 
can be decomposed as 



H - 22 ^^^^ + zJ ^i 



k.Kn 



i<n 



where H^^i is a self-adjoint operator acting on the joint 
Hilbert space of particle k and / and Hj is the free Hamil- 
tonian of particle j (without loss of generality, we can 
drop the second sum, by reckoning it to the first part). 
Pair-interactions Hk^i are the infinitesimal versions of 
two-qubit gates: every unitary of the form e'^'-''* for 
t € R is a two-qubit gate. On the one hand, this seems 
to be an important justification for two-qubit gates, since 
it refers to the form of the fundamental forces of nature, 
on the other hand, the argument is not really correct: In 
general there is no obvious correspondence between the 
time evolution 



and any finite sequence of two-qubit gates. However, 
there is an obvious simulation by two-qubit gates in an 
approximative sense given by the well-known Trotter for- 
mula: 



lim (Hfcje*^'=''*/")" 
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This example shows, that the simulation of the time 
evolution caused by a time-independent pair-interaction 
Hamiltonian might require an infinite number of two 
qubit gates. Hence it seems to suggest that a defini- 
tion of complexity based on two-qubit gates does not take 
into account the most natural form of dynamics of man; 
particle quantum physics. However, it has been shown 
that the number m of gates required to simulate the right 
hand side of (0) up to an error e is only growing with i^/e. 
Despite the fact that infinite accuracy requires an infinite 
number of gates, the time for implementing the growing 
number of unitaries does not tend to infinity if one as- 
sumes that the implementation of e^Hk,it/m j-gq^jj-eg t:he 
time 0{t/m) (see [||). Taking this assumption, we will 
show that the running time of a discrete quantum algo- 
rithm for simulating the dynamics given by (|l|) is deter- 
mined by the chromatic index of the interaction graph. 
This turns out to be true even for time-dependent pair- 
interactions. One might reformulate this result by saying 
that the running time depends on 'the complexity of the 
interaction'. The relevance of this result is twofold: In 
case realizations of future quantum computers are based 
on two-qubit gates, it gives exact statements about the 
complexity of simulating non-autonomous quantum dy- 
namics resulting from pair-interaction Hamiltonians. To 
our knowledge, this is the first exact analogue between 
complexity measures of discrete and continuous quantum 
algorithms. Secondly, our result is relevant for the sim- 
ulation of arbitrary pair-interaction dynamics by a dy- 
namic with restricted interaction graph: We derive state- 
ments about the efficiency for simulating time evolutions 
based on pair-interactions with high chromatic index by 
other evolutions with interactions of lower index. 
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II. DISCRETE AND CONTINUOUS QUANTUM 
ALGORITHMS 

In our model of discrete quantum computers, two-qubit 
gates acting on disjoint pairs of qubits can be imple- 
mented simultaneously. We define: 

Definition 1 A discrete quantum algorithm A of 

depth k is a sequence of k steps {Ai,...,Ak} where 
every step consists of a set of two qubit gates {uji}j,i 
acting on disjoint pairs (j, I) of qubits. Every step i 
defines a unitary operator Vi by taking the product of 
all corresponding unitaries in any order. The product 
u := Iii<k'Vi is the 'unitary operator implemented by A'. 

The following quantity measures the deviation of a uni- 
tary operator from the identity: 

Definition 2 The angle of an arbitrary unitary opera- 
tor u is the smallest possible norrnn \\a\\ of a selfadjoint 
operator a which satisfies e*"^ — u. It coincides with the 
time required for the implementation of u if the norm of 
the used Hamiltonian is 1. 

This term allows us to formulate a modification of the 
term 'depth' which will later turn out to be decisive in 
connecting complexity measures of discrete and continu- 
ous algorithms: 

Definition 3 Let ai be the maximal angle of the uni- 
taries performed in step i. Then the weighted depth is 
defined to be the sum ^^ ai. 

Assuming that the implementation time of a unitary 
is proportional to its angle, the weighted depth is the 
running time of the algorithm. Since this coincidence is 
based on a possibly unrealistic assumption we will prefer 
the term 'weighted depth' in formulating exact mathe- 
matical statements. 

Now we want to formalize the notion of quantum al- 
gorithms based on time dependent Hamiltonians. Such 
continuous algorithms have already been considered in 
the hteratureQ (e.g. see |j-||]), but in our approach the 
Hamiltonians are explicitly restricted to pair-interactions 
in many-particle systems: 

Definition 4 A continuous quantum algorithm A 

of running time T is a piecewise Lipschitz continuous 
function t h- > H{t) from the interval [0, T] into the set of 
those self- adjoint operators acting on 7i„ which can be 
decomposed into 



where every Hk,i is a self-adjoint operator acting on the 
qubits k and I. Here Hk^i = Hi^k imd the sum runs over 
the unordered pairs. We say 'A implements u' if u — ut 
and {ut)teT is the solution of the non- autonomous dif- 
ferential equation -^Ut — —iH{t) with uq = 1. 

As already mentioned above, the weighted depth of a 
discrete algorithm for simulating the evolution (|l|) is fi- 
nite. More precisely, it is at most mt if m is the number 
of pairs (fc, I) with the property that H^j ^ 0. But, if 
one has additional information about the interactions, 
one can make further statements about possible paral- 
lelization. This is illustrated by the following examples: 
In case all the pairs (fc, Z) with H^^i ^ are disjoint, we 
can perform all the unitaries e*^*"'*'™ simultaneously and 
get the running time t. In case of the nearest-neighbor 
interaction in a one dimensional spin chain, we have only 
pairs of the form (fc, k + l). Hence we can perform all the 
transformations 

gi-H'2fc,2fc + l*/»n 

simultaneously and the operations 

ii?2fc+l,2fc + 2*/m 

as well. Here, parallelization allows to decrease the run- 
ning time down to 2i. 

We show that the graph theoretical concept of chro- 
matic index offers the appropriate terminology for deter- 
mining the degree of possible parallelization. 



H{t) =Y,HuAt), 



k,l 



Definition 5 |Jq/ The chromatic index of a graph 
with vertices {1, ... ,71} is the least number of colors re- 
quired for coloring the edges in such a way that there are 
no edges with the same color having a common vertex. 

For every time t during the running time of the contin- 
uous algorithm we define a family of undirected graphs: 

Definition 6 For every non-negative real number r and 
every time t S [0,T] we define the interaction graph 
Gr{t) as follows: Take the qubits {l,...,n} as vertices 
and let the edges be all the pairs (fc, I) with the property 

\\HkAm>r. 

The following quantity turns out to be decisive for the 
degree of possible parallelization of the discrete simula- 
tion: 



^Here and in the following ||.|| denotes the operator norm given by ||at| := max^, ||aa::|| where x runs over the unit vectors of 
the corresponding Hilbert space. 
■^In W time-dependent Hamiltonian algorithms are called an 'analog analogue of a digital quantum computer'. 



Definition 7 Let nr{t) he the chromatic index of Gr{t). 
For any time t € [0, T] of a continuous algorithm, we de- 
fine the weighted chromatic index W{t) as 



W{t) 



nr{t)dr. 



Furthermore we will need the time integral of W: 

Definition 8 The integrated chromatic index of the 

quantum algorithm is defined to be 

/ := / / nr{t)drdt, 
Jo Jo 

where T is the running time of the continuous algorithm. 

Actually, the terminology 'integrated weighted chro- 
matic index' would be appropriate. We preferred the 
short terminology. 



III. TRANSLATION BETWEEN DISCRETE AND 
CONTINUOUS ALGORITHMS 

The following theorem suggests that the simulation 
of an interaction with high integrated chromatic in- 
dex generically requires a discrete algorithm with high 
weighted depth. More precisely, it shows that the com- 
plexity of a unitary transformation can equivalently be 
defined as the infimum over the values of the weighted 
depth of all possible discrete implementations or the in- 
fimum over the values of the integrated chromatic index 
of all possible continuous implementations: 

Theorem 1 

a) Every arbitrary unitary operator u acting on n qubits 
which can be implemented by a discrete quantum algo- 
rithm with weighted depth a can also be implemented by a 
continuous quantum algorithm with integrated chromatic 
index a. 

b) If there is a continuous algorithm with integrated chro- 
matic index a implementing u then there is a sequence of 
discrete algorithms (Ai;)^^^ implementing the unitaries 
Uk such that limfe_,oo Uk = u and the corresponding val- 
ues of the weighted depth converge to a. 

Proof: Statement a) is almost trivial: Let A be a 
discrete quantum algorithm of depth k with the maximal 
angles {aj)j<k- Let Sj be the set of qubit pairs which are 
addressed in the step j and {up)p(zs^ the corresponding 
set of unitary transformations. Define a continuous algo- 
rithm with running time fc as follows: for t G (i~l,j)weH 
define the constant Hamiltonian H{t) :— X]ne5 ~*lnup. 



Since the angle of a unitary operator u is the norm of 
Inu, we have aj — maxpg^. || lnup||. Hence, for every t 
in the interval (j — 1, j) we have nr{t) = for r > aj and 
nr{t) = 1 otherwise. Obviously we have 



nr{t)drdt ~ 2_, ^j- 

3 



For statement b) take a partition of [0, T] into small 
intervals on which H is continuous. Let (i,t + e) be one 
of those intervals. Let 



ri < 



< r„ 



be the positive numbers {\\Hij{t + e/2)||)ij<„ in its 
canonical ordering. Furthermore define ro := 0. Obvi- 
ously, the function r i-^ nr{t-\- e/2) takes constant values 
Uj on every interval ]rj-i,rj] for 1 < j < to. Let 

Af 1 U M2 U . . . U M„^. 

be a partition of the set of edges of Gr{t + e/2) given 
by an allowed coloring (in the sense of Definition 0) of 
the edges corresponding to the chromatic index hj. For 
every j < to we proceed as follows: For every Mp define 
a step of a discrete algorithm by the set of unitaries 



|e(-H- 



fc,,(t+e/2))x(rj-rj_i)/||Hfc,i(t+£/2)|| 



{k,l)£Mp 



Since p runs from 1 to rij we obtain Uj steps for the 
discrete algorithm. Every step has at most the angle 
e{rj — rj_i). The weighted depth of these Uj steps is 
_i). Doing this for every j < to, we 



at most nje(rj 



obtain a sequence of J2j<m'"'j steps which substitutes 
the continuous algorithm on the interval (i, i -I- e) up to 
an error in the order of e^. The weighted depth of this 
sequence is at most 



Yl '^J^'-^'J 



rj-i 



1<J?<?71 



nrit + e/2)dr. (3) 



Without loss of generality we will assume the weighted 
depth to be equal to the terms in equation (^) since we 
can blow up the algorithm by transformations which can- 
cel each other out. If the total discrete algorithm is de- 
fined by combining the sequences for every interval of 
the form (s, s -|- e) we get a simulation of the continuous 
algorithm with an error in the order of e. Furthermore 
the total weighted depth converges to L L nr{t)dt for 
e^O. □ 

One can consider discrete quantum algorithms as spe- 
cial cases of continuous quantum algorithms in the sense 



^Here we denote an interval by (.) if we do not care about whether it is open or closed. 



of the proof of part a) of Theorem |l|. Then the dis- 
crete quantum computer is obtained by restricting the 
Hamiltonians H{t) to those with chromatic index one 
and norm one for each non- vanishing component Hk^i(t). 
In this sense, part a) of Theorem n^ is a special case of 
the general principle that continuous time algorithms can 
be simulated by other continuous time algorithms with 
lower chromatic indices and the same strength of the pair- 
interactions where the running time is increased by the 
quotient of the chromatic indices. We sketch this simple 
observation in a not too formal way: 
At time t we have the interaction 



H{t) 






where (j, I) runs over the edges E of the graph Go{t). For 
m < no(t) take k such that km > nQ(t). Take a parti- 
tion of the edges of Go{t) into k subsets such that the 
subgraphs with the set of edges {Ei)i<k have chromatic 
indices less or equal to m. Substitute the small time in- 
terval (t, i + e) of the original algorithm by the k time 
intervals (t + (i — l)e, t + ie) with 1 < j < fc in which the 
Hamiltonian J^a 06_e ^j,ii^ + ^/2) is switched on. The 
difference of the unitary operators implemented by the 
substituted algorithm and the original one is of second 
order in e for every subinterval. Hence the total error 
tends to zero for e ^ 0. 

This observation justifies in some sense our point of 
view that one might think of the chromatic index as the 
complexity of the interaction. This raises the question 
whether bounds can be given on the required integrated 
chromatic index necessary for preparing certain entan- 
gled quantum states from an initial product state. 



IV. COMPLEXITY BOUNDS FOR HIGH 

ENTANGLEMENT 

In ]ll[ | we gave lower bounds on the depth required 
for preparing certain states with large-scale entanglement 
(some of those bounds are easy conclusions from [Q ) for 
large qubit numbers n. More specifically, one can show 
the following: Let {aj)j<n be a family of arbitrary self- 
adjoint operators where every Uj acts on qubit j and has 
operator norm 1. For product states, the variance of the 
observable 



'■= E °J 



(4) 



j<n 



grows with 0{n). Values of the variance in the order of 
n^ indicate large-scale entanglement for pure states [Oj . 
In a state obtained by a discrete algorithm with depth 
k we could show |11| the variance to be less or equal to 
n 4*^/2. Hence the emergence of large-scale entanglement 
in the sense described above requires a depth in the order 
of logn. 



It is natural to ask whether similar bounds can be 
shown for the weighted depth and the integrated chro- 
matic index. However, below we present a proof for 
such a bound. In contrast to the discrete definition of 
depth the bound is asymptotically independent on n in 
the limit n — > oo. Whether this is a lack of our estima- 
tions or whether it is a hint for a fundamental difference 
between complexity measures of continuous and discrete 
algorithms is unclear. We have: 

Theorem 2 Let t h^ H{t) be a continuous quantum al- 
gorithm with integrated chromatic index a implementing 
the unitary u. Let the quantum computer start in the 
product state p®". Let a :— u{p®^^u^ be the state ob- 
tained by the algorithm. Then for the variance of any 
observable a of the form M) in the state a we have 



V^a) < 



(1 - 2ay 



with 



Va{a) := tr[aa^) - {tr[aa)f 

= tr{p'^'\u^au)iu^au)) - tr{p'^''{u^au)f. 

Proof: Due to Theorem |l| we can restrict the proof 
to the case that H is of the same form as that one con- 
structed in the proof of Theorem |l| a), i.e., H{t) is for ev- 
ery t a sum of pair-interactions acting on disjoint pairs. 
Furthermore we will assume without loss of generality 
that the norm of every non-vanishing pair-interaction is 
one. Then a coincides with the running time T. We write 
H{t) = J2ei£E(t) ^ei^) where E{t) is the set of edges of 
the interaction graph Go(i). By solving the differential 
equation 



dt 



a{t) = i[H{t),a{t)] 



with a(0) = a we obtain u^au = a{T). Explicitly, u'^au 
is given by a Dyson series |14| : 

u^au = E *' fiH{tk), [. . . , [ff(ii), a] . . .]]dt\ 

feGN •' 

where the integration is carried out over the simplex 

<ti < ... <tk <T. 

Let X :— (pk, . . . ,pi,r) he a, k + I tuple consisting of 
k unordered pairs pj (which are considered as subsets of 
{1, . . . , n}) and one element r G {1, . . . , n}. Define 

Bx{tk, ■ ■ ■ ,ti) '■= 

[HpSk), [Hp,_,{tk-i), [..., [Hp,{ti),ar] . . .]]. 

Furthermore, for an arbitrary pair of observables c, d we 
introduce the covariance C{c, d) in the state p®" as 



C{c,d) := tr{p®''cd) - tr{p'^'^c)tr{p''^''d). 



Hence we get 



Then the variance of the observable a in the state a 
can be written as: 

k,l 



with 



Fky.= I lY,C{Bx{tk,.-.M)BY{ 



Si,.. . ,si))dt^ds , 



where the sum runs over aU {k + l)-tuples X of the form 
X := (cfc, . . . ,ei,r) with Cj G E(tj)\/j and r G {1, . . . , n} 
and (Z + l)-tuples Y which are defined analogously by us- 
ing the sets of edges E{sj) instead of E{tj) (here k runs 
over all non- negative integers) . In the following, we shall 
sometimes consider an X of this form canonically as a 
subset of {1, . . . , n} and write j G X iff j agrees with the 
rightmost element of X or with one of the vertices of at 
least one e^ for 1 < j < fc. 

Now let the times tk, . . . ,ti and s;, . . . , si be fixed. We 
define Xk as the set of (fc + l)-tuples (efc, . . . , ei, r) with 
Cj G E{tj) and the property that every Cj has a vertex 
which is an element of (ej_i, . . . , ei, r). It is easy to see 
that Bxitk,. . . ,ti) ^ implies X e Xk- The set X^ 
can be constructed iteratively: Set Xq := {(1), . . . , (n)}. 
Then we have (ek, X) G Xk if and only if X G Xk-i and 
Cfc G E{tk) and at least one of the vertices of et is an el- 
ement of X (here (e^, X) is an informal notation for the 
(fc-l-l)-tuple obtained by appending Ck to the fc-tuple X). 
In an analogous way we define the set yi by referring to 
the sets E{sj) instead of E{tj). 

The covariance 



C{Bxitk,...,h)BYisi,...,. si)) 



(5) 



can only be nonzero ii X G Xk and Y E yi 
have nonempty intersection (considered as subsets of 
{1, . . . , n}). We show that there are at most 
n{k + 1)!(Z -|- 1)! pairs X, Y with this property by proving 
that for every j G {1, . . . , n} there are at most (fc + 1)! 
sets X £ Xk with j £ X. This can be seen by induc- 
tion over fc: For fc = the statement is obvious. Assume 
the statement to be true for Xk-i- For every X G Xk-i 
with j d X we have at most fc -I- 1 elements of the form 
{ek,X) G Xk since for each vertex v E X there can be 
at most one ek G E{tk) such that Ck has w as a ver- 
tex (remember that the corresponding graph GQ{tk) has 
chromatic index 1). If j ^ X the statement j G {ek,X) 
can only be true if j G ek G E(tk). There is at most 
one Cfc G E{tk) with vertex j. Let ek be the pair {j, m}. 
By assumption there are at most fc! tuples X G Xk-i 
such that m G X. Hence there are at most (fc -f- l)fc! sets 
(efc,X) €Xk with j G (efc,X). 

Since the norm of every i?e(i) is zero or one by as- 
sumption, the norm of Bx (t/t , . . . , ti ) can never exceed 
2fc_ 



J2x£Xk.Yeyi C{Bx{tk,-. .,ti),BY{si, 



,si)) 



Since the integration is carried out over a simplex of size 
TV /{kill) we obtain: 



K(a) <n}^{2T)^{k + l){2Ty{l + l) 

k,l 



n(V(2r)^-(fc + i))' 



This power series converges for 2T < 1. Since 

^^(fc + l)p^ = 1/(1 -p)^ for every - 1 < p < 1 



fceN 



we have: 



K(a) < 



(1 - 2Ty 



Theorem ^ shows that for huge n the variance of the 
measurement results of any mean-field observable can 
only grow in the order of n^ if the integrated chromatic 
index a is at least 1/2. Hence a > 1/2 is a necessary 
condition for the particular class of large-scale entangled 



states considered in |l3 



V. CONCLUSIONS 

We have shown that the running time of a discrete 
algorithm simulating the evolution defined by a time- 
dependent Schrodinger-equation with pair-interaction 
Hamiltonians is given by the time integral over the chro- 
matic index of the interaction graph. This result sug- 
gests to take this time integral as a complexity measure 
for continuous time algorithms and it seems natural to 
ask for the complexity for the preparation of highly en- 
tangled states in many-particle systems. We proved such 
a bound. At the moment, we cannot decide whether it 
is tight or not. Given a quantum system with the prop- 
erty that the weighted chromatic index W{t) of the rele- 
vant interaction Hamiltonian satisfies W{t) < ni for ev- 
ery time t we conclude that the system requires at least 
the time l/{2m) in order to produce large-scale entan- 
glement in the sense explained above. If we assume a 
non-vanishing probability of depolarizing errors for each 
single qubit, the error probability has to decrease with 
0(l/n) if large-scale entanglement should be maintained 
||l3| . This strongly suggests that the chromatic index or 
the strength of the interaction has to increase in the order 
of n. 

For nearest-neighbor interactions in solid state physics 
it is not useful to apply our estimations, since there are 



considerably tighter bounds (see [|5[ paragraph 6.2.1) 
But for mean- field interactions Jl^], where the weighted 
chromatic index W{t) is already determined by the 
strength of the pair-interactions, we cannot see any ob- 



vious tighter bounds. Hence physical systems like solid 
states with long-range interactions will present a typical 
application of our results. 
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